Ensemble dependence in the Random transverse-field Ising chain. 
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In a disordered system one can either consider a microcanonical ensemble, where there is a precise 
constraint on the random variables, or a canonical ensemble where the variables are chosen according 
to a distribution without constraints. We address the question as to whether critical exponents in 
these two cases can differ through a detailed study of the random transverse-field Ising chain. We 
find that the exponents are the same in both ensembles, though some critical amplitudes vanish 
in the microcanonical ensemble for correlations which span the whole system and are particularly 
sensitive to the constraint. This can appear as a different exponent. We expect that this apparent 
dependence of exponents on ensemble is related to the integrability of the model, and would not 
occur in non-integrable models. 

PACS numbers: 05.60.-k, 72.10.Bg, 73.63. Nm, 05.40.-a 



I. INTRODUCTION 

In the study of the critical behavior of disordered sys- 
tems, it is usual to pick the random variables from some 
distribution. This allows sample-to-sample fluctuations 
in the sum of the interactions (e.g. nearest neighbor) of 
order y/~N. We will call this the canonical ensemble of 
disorder, by analogy with the canonical ensemble of sta- 
tistical mechanics which allows fluctuations in the energy. 
It is sometimes of interest to complete this analogy and 
define a microcanonical ensemble of the disorder in which 
there is a strict constraint, for example by fixing exactly 
the sum of the (e.g. nearest neighbor) interactions in 
each sample. Our experience from conventional statisti- 
cal mechanics tells us that in the thermodynamic limit 
the choice of ensembles does not matter but it is not very 
clear that this is also true for random systems, especially 
for the case of quantum phase transitions. 

In this paper we will study the simplest disordered 
model with a quantum phase transition, the random 
transverse-field Ising chain (RTFIC)^ with the Hamil- 
tonian: 
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where Jj > and hi > are random variables cho- 
sen from distributions p(h) and 7r(J) with averages 
[ln/i] a v, [In J] av and variances var(ln/i), var(lnJ). We 
use free boundary conditions, so the sum for the J, stops 
at L- 1. 

Let us define the two ensembles, microcanonical and 
canonical, precisely for this model. For the canonical 
ensemble the hi and Jj are chosen randomly. A parame- 
ter which characterizes the deviation from criticality is 8 
where 



[In hi 



[InJL 



(2) 



var(ln h) + var(ln J) ' 
For the microcanonical ensemble we constrain the hi and 



Ji such that the parameter 



1 Z)f=l ln hi ~ ln J i ~ Pg h h 

L — 1 var(ln h) + var(ln J) 
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is set to a prescribed value for each sample. The last 
term in the numerator, which is not necessary to get the 
asymptotic behavior, corrects for there being one more 
hi than Jj with free boundary conditions. It ensures 
that 8 = [(5] av even for a finite system. In the canonical 
ensemble, the fluctuations in 8 from sample to sample 
are 0(1/ VT). It is known that a phase transition occurs 
in this model at 8 = 0. All our numerical results in the 
paper will be at the critical point. 

Pazmandi et al£ have argued that it is precisely the 
0(1/ y/L) fluctuations which lead to the bound on the 



finite size correlation length exponen 
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> 2/d. Fur- 



ther, they argue that exponents in the microcanonical 
ensemble need not satisfy this bound. In addition, Igloi 
and Riegeri claim that the exponents of the RTFIC can 
depend on the choice of ensemble. In this paper our 
main goal is to investigate these claims through a de- 
tailed study of the zero-temperature critical properties 
of the RTFIC. 

Using the Jordan- Wigner transformation, the Hamil- 
tonian of the RTFIC can be mapped to a free-fermion 
problem and this mapping is particularly useful in the 
context of numerical computations. It can be shown that 
various physical quantities can be expressed in a straight- 
forward way in terms of the eigenvalues and eigenstates of 
the free Hamiltonian, which are easy to evaluate numer- 
ically. These have been discussed by several authors in 
earlier papers^ and we will not repeat the derivations 
here but will use those results in our numerical calcu- 
lations. In this paper we look at the surface and bulk 
magnetizations, the end-to-end correlation function and 
the energy gap. For the surface magnetization, the free 
fcrmion method leads to a simple form and it is possible 
to obtain some detailed results analytically. We first dis- 
cuss this in Sec. Inland then present numerical results for 
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various other quantities in Sees. Illll and II VI Our conclu- 
sions are summarized in Sec. [V] 



II. SURFACE MAGNETIZATION 

The simplest quantity to calculate is the surface mag- 
netization, m s , which is defined with free boundary con- 
ditions in which we fix a z at one end, i = L say, to be +1. 
The surface magnetization is then the expectation value 
of <7 Z at the other end (i = 1). This is equivalent to delet- 
ing the transverse field on site L, so a z L commutes with 
the Hamiltonian and the ground state is exactly doubly 
degenerate, and calculating the expectation value of a\ 
in the ground state with ct£ = 1 . Let us denote this state 
by |0) and so 



This has a simple formic, namely: 
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Igloi and Riegeri used this to numerically compute the 
distributions of B = — log(m s ), P C {B) and P mc (B), in 
the two ensembles. 

However it is also possible to obtain the distribution 
functions analytically 10 for L — > oo and we rederive those 
results here. First consider the canonical case. Let xi = 
(hi/Ji) 2 . Then from Eq. © we get 

B L = ^ln[l+x 1 (l+x 2 {l+x 3 (l + ...{l+x L - 1 )...)))} 



= -ln[l+xie 
« ^ ma; i + b l-i 



(6) 



The above approximation is good most of the time since 
£>l is expected to be order L 1 / 2 . For small Bl we notice 
that the increment in Bl is always positive and so Bl 
can never become negative. This and Eq. © means that 
Bl can be effectively described by a biased random walk 
(in which L is the time variable) with a reflecting wall at 
the origin. It is convenient to introduce a scaled length 
variable 



£ = L 



var(ln h) + var(ln J) 



(7) 



in terms of which the probability distribution P C (B,L) 
can be written as 



P C (B,L) = P(BJ). 



(8) 



Then in the continuum limit, it is easy to see that P(B, £) 
satisfies the following equation: 
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where 5 is given by Eq. J5J • The reflecting boundary con- 
dition is imposed by requiring the current at the origin 
B = to be zero, thus [ff - 2SP] B =o = 0. This prob- 
lem is mathematically equivalent to Brownian motion in 
a gravitational field and its solution, for identical bound- 
ary conditions, is discussed in Ref. [TJ With the initial 
condition P(B, L = 0) = 5(B) we find that the solution 
of the above equation isi -: 



P C (B,L) = P(B,£) 



(10) 
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where erfc is the complementary error function. 

For the microcanonical ensemble the distribution 
P mc (B, L) can be found using the result that it is related 
to P C (B, L) through the general transformation Eq. I|A6() . 
One then get si -: 



P""'(D.L) = 20(B) (-){B - 25l ) [j-5)< 



~2BS 



(11) 

Note that P C (B, L) is a function of the two scaling 
variables b = B/i 1 ' 2 and U 1 ' 2 , and similarly P mc (B, L) 
is a function of 5l x l 2 as well as b. According to finite size 
scaling, the scaling variable associated with the deviation 
from criticality (5 or i5 here) should be proportional to 
L 1 /". Hence Eqs. <jTU|) and (fTT|) show that the true cor- 
relation length, as determined from finite size scaling, is 
v = 2. 

Now that we have the complete distributions for B = 
— ln(m s ) we can calculate the mean surface magnetiza- 
tion. Even though we find that [— ln(m s )] av ~ L 1 / 2 in 
both ensembles, the behaviour of the mean of m s (rather 
than its log) is quite different in the two ensembles. This 
is because m s — e~ B and so the main contribution to 
m s comes from the behaviour of P(B) at small B (i.e. 
from rare samples with m s sa 1). For large L we find the 
following asymptotic forms for the mean magnetization. 
In the canonical case: 



(ni) 1 / 2 



5erfc[5£ 1/2 } 



(12) 



while in the microcanonical ensemble we get: 



= -r-25 



(6<0) 



-281 



(25+-) (5>0). 



(13) 



~ exp(-£/£ c ) 
= 2 deduced ear- 



For 51 1 ! 2 > 1, Eq. JHJ gives [m 
where £ c = \/5 2 , in agreement with v 
lier. However, for the microcanonical distribution with 
51 > 1, we have [m s ]Z ~ exp(-£/£ mc ) with £ mc = 1/5. 
This looks like an apparent correlation length exponent 
of 1 rather than 2. However, it is worth investigating 
the origin of this discrepancy between the apparent ex- 
ponents in the two ensembles. For both ensembles the 
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scaling variable is 5 (or 5) L 1 / 2 . In the canonical ensem- 
ble, the distribution in Eq. (|l(Jf> has a constant weight at 
B = 0, which leads to the expected £ c = 1/5 2 . However, 
for the microcanonical ensemble, there is a "hole" in the 
distribution for B < 25£. The average is dominated by 
the part of the distribution with the smallest B, so this 
difference in the distributions for small B accounts for 
the difference in the behavior of the average. Since the 
weight of the distribution for the microcanonical ensem- 
ble vanishes at small B, we argue that the amplitude of 
the expected 1/5 2 divergence of the correlation length for 
[m s ]™ c is zero, and that the resulting behavior, ~ 1/5, 
is really a correction to scaling. In the rest of this paper 
we shall reinforce the conclusion that v = 2 in both en- 
sembles but with the leading amplitude vanishing, in the 
microcanonical case, for certain quantities which are par- 
ticularly sensitive to the microcanonical constraint. This 
point of view is different from that of Igloi and Riegep£ 
who argue that v is different for the two ensembles. 

For L — > oo and 5 (or 5) < 0, we get [m s ] av = — 25 in 
both ensembles, giving a magnetic exponent (3=1. At 
the critical point we find that the mean magnetization 
decays with system size as: 
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FIG. 1: Mean of the bulk magnetization, as defined Eq. I16L 
for different system sizes, evaluated in the two ensembles. The 
expected slope from finite size scaling is j3/v ~ 0.191. 
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From finite size scaling we expect a decay ~ L~PI V where 
v is the correlation length exponent. While this might 
suggest that v = 2 in the canonical ensemble and v = 1 
in the microcanonical case, we feel, as discussed above, 
that a more consistent picture is that the amplitude of the 
leading divergence of the correlation length appropriate 
to [to s ]™ c vanishes for the microcanonical ensemble and 
that the true exponent is v = 2 in both cases. 



middle of the chain 



m = (6K /2 |o). 



(16) 



Using the free fermion method this can be expressed as 
the determinant of a matrix whose elements are expressed 
in terms of eigenstates of a quadratic Hamiltonian. We 
evaluate this numerically and compute both the mean of 
the bulk magnetization and also its distribution in the 
two different ensembles. Here we will only examine the 
data at the critical point. 

Since the spin at one end is fixed, there are equal num- 
bers of hi and Jj, so the definition of 5 in Eq. © is slightly 
modified, which leads to the condition 



III. BULK MAGNETIZATION 

In the previous section we saw that the correlation 
length exponent for the surface magnetization seems, at 
first glance, to be different in the canonical and micro- 
canonical ensembles, but we argued that the correct in- 
terpretation is that the exponents are the same, v = 2, 
but the amplitude of the expected divergence of certain 
quantities is actually zero for the microcanonical ensem- 
ble. In this section we strengthen this argument by in- 
vestigating the magnetization in the bulk of the sample, 
when a spin at the end is fixed. We see the same ex- 
ponent v = 2 in both ensembles, clearly indicating that 
a correlation length with a v = 2 divergence does exist 
for the microcanonical ensemble. Its absence in the sur- 
face magnetization presumably indicates that the leading 
amplitude vanishes for this quantity. 

We again consider an open chain with the spin at one 
end fixed to o~ z L = 1 and look at the magnetization at the 
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(17) 



for criticality (i.e. 5 = 0) in the microcanonical ensemble. 
We set J = 1 and allow h to take values 2 and 1/2. In 
the canonical case each hi takes one of its two values with 
equal probability. In the microcanonical case exactly half 
of the sites, chosen at random, are assigned h = 2 and 
the other half are assigned h = 1/2, which clearly satisfies 
Eq. QZJ. 

The numerical results for the decay of the mean mag- 
netization with system size are shown in Fig. ^ The 
mean is seen to behave similarly in both ensembles and 
the system-size decay is consistent with the form ex- 
pected from finite-size scaling [m] av ~ L~PI V with j3 = 
(3 - 5 1/2 )/2 and v = 2, so that /3/v ~ 0.191. 

We now look at the distribution of to. We use the 
variable b = — \n(m)/ L 1 / 2 since this has good scaling 
properties. The details of the distributions of to, shown 
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FIG. 2: The distribution of the bulk-magnetization for dif- 
ferent system sizes for the canonical ensemble. 

in Figs. [21 and^l are different in the two ensembles; in 
particular the microcanonical distribution falls off faster 
at large argument. However, as can be seen in Fig. 01 the 
behaviour at small values of the argument is the same 
in both ensembles, which leads to the same asymptotic 
behaviour for [m] av . 

Thus, unlike the surface magnetization, the bulk mag- 
netization shows the same critical behaviour in both en- 
sembles. In particular, the results of this section indicate 
that there is a correlation length which diverges with ex- 
ponent v — 2 in the microcanonical ensemble. It is not 
seen in the surface magnetization m s , for which the cor- 
relation length diverges less strongly with an exponent 
v = 1, but this must simply indicate that the amplitude 
of the v = 2 divergence vanishes for m s . 



IV. END-TO-END CORRELATIONS AND GAPS 

In this section we investigate numerically the energy 
gap A and the end-to-end correlation function 



c 1L = (okxio) 



(18) 



in the canonical and microcanonical ensembles to com- 
pare the results of the two ensembles with each other and 
with analytical results^ for the canonical ensemble. 

We take the following rectangular distribution for the 
bonds and fields at the critical point 



tt(J) 



p{h) = 



1 for < J < 1 

otherwise 

1 for < h < 1 
otherwise. 



(19) 



which gives 



FIG. 3: The distribution of the bulk-magnetization for dif- 
ferent system sizes for the microcanonical ensemble. 
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FIG. 4: Comparison of the distributions of the bulk- 
magnetization in the two ensembles for a chain of length 
L = 201. On this log-log plot they appear to have the same 
slope at small values of the argument. However, at large ar- 
gument, the microcanonical distribution falls off faster. 



[lnJ] av = -1, var(lnJ) = l. (20) 
From this it follows that 

e = L, (21) 

where I is defined in Eq. Q. We use free boundary 
conditions without constraining either of the end spins. 
From Eqs. © and 1)20(1 the condition for criticality in the 
microcanonical ensemble is 
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FIG. 5: Plot of the mean values of the end-to-end correlation 
function Ci,l and the energy gap A for both the canonical and 
microcanonical ensembles. In the canonical case, the fits are 
close to those predicted analytically and given in Eqs. J33 
and 12 11 . For microcanonical ensemble, the lines are just 
guides to the eye. 



We initially generate the hi and Ji in an unconstrained 
way, as for the canonical ensemble, but then rescale hi 
by an appropriate factor so that the above condition is 
satisfied. 

The distributions of C\l and A were studied earlier by 
Fisher and Young 2 and we summarize some of their main 
results for L — > oo: 



[Cil]\ 
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ln(A)-ln(C 1L ) = ^ln^O-^MJ/; 



(24) 



(25) 



, ver- 
2 , for 



FIG. 6: Plot of the mean correlation function [Ci,l 
sus system size compared with [mfrofjav an d ([Tif]a 
the canonical ensemble. The ([mf] av ) 2 data falls as l/L as 
expected from Eq. 114H . 
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In Fig. we compare the system size dependence of 
the average correlation function and the gap in the two 
ensembles. For the canonical case, the data agree well 
with the analytic predictions in Eqs. and 124(1 . as 
was also found in Ref. 0- However, if we fit the data for 
the gap to Eq. I|24|) adjusting only the overall amplitude 
the x 2 is 150 which is very high. Hence there must be 
systematic corrections to Eq. (|24|) which are larger for 
the sizes studied than the, very small, statistical errors. 

For the microcanonical ensemble, the data for both 
[Cil]w° an d [A]av C in Fig. appear to decay as stretched 
exponentials, and we will discuss fits to this data below. 

Interestingly, the value of [Ci rj a v is found to be close to 
[mfmfj av , where raf (m s L ) is the surface magnetization 
at site 1 (L) with the spin at site L (1) fixed. This can be 



FIG. 7: Plot of the mean correlation function [Ci,L]av C versus 
system size compared with [mf mfj™ c and ([mf]™ c ) 2 , for the 
microcanonical ensemble. The data for ([mj]™) 2 decays as 
l/L 2 as expected from Eq. 1151 . 



seen in Figs. and [7| where we plot both these quantities, 
as well as [mf] 2 v , for different system sizes. Especially 
in the canonical case, [CiL] av and [mlm s L ] a _ v are almost 
indistinguishable. Note that for the microcanonical, but 
not the canonical, ensemble [m^] 2 v is much greater than 



l LJav 



The reason for this is that 



l lJav 



is dominated 



by a few rare samples where the bonds are bigger than the 
fields at the free end (and hence, for the microcanonical 
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FIG. 8: Plot of [m{m 3 L ]w C against L for different system 
sizes, for the microcanonical ensemble. A best fit to the form 
aL p e- bLl/3 with p = 0.44 (x 2 = 3.9) is shown. The inset 
plots x 2 against p. 
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FIG. 9: Distribution of [In (A) - ln(Ci L ) + 1]/(2L) 1/2 for the 
canonical case. The analytic form, deduced from Eq. 1251 . is 
a Gaussian with variance unity. This is shown by the solid 
line. 



ensemble, must be less than the fields at the other end 
because of the constraint). Hence for these samples m s L 
is smaller than typical in the microcanonical ensemble. 

Since [Cix] av and [mfm s L ] a _ v behave similarly we can 
obtain a better estimate for the decay law in the micro- 
canonical case. This is because mlm B L can be obtained 
directly from Eq. JSJ and thus be accurately computed 
numerically for bigger systems. We assume the same 
form as in the exact result for the gap in the canoni- 
cal case, Eq. (|24[) . i.e. aL p exp(— bL^) : and take /j, = 1/3 
the same value as in Eq. I|24|) . The data shown in Fig.[S]is 
fitted to aL p exp(— 6L 1/3 ) by varying a and b for several 
(fixed) values of p. The minimum \ 2 of 3.9, which is quite 
acceptable for three degrees of freedom, is obtained for 
p ~ 0.44. It we assume that p = 1/6, as for the canonical 
case, then \ 2 = 1380 which is extremely high. However, 
we noted for the canonical case, that there appear to be 
corrections to the scaling form in Eq. I|24(l . Hence we 
cannot rule out the possibility that p = 1/6 also for the 
microcanonical ensemble. 

From Figs. and it seems plausible that 
[mfmfj™ , [Cil]™ c and [A]™ c all vary in the same way 
in the microcanonical ensemble. If this is so then the 
data for the gap is consistent with the stretched expo- 
nential form aL 1 / 6 exp(— 6i 1/3 ), for both canonical and 
microcanonical ensembles, though there are some sys- 
tematic corrections to this for the range of sizes that can 
be studied. This is known to be exact for the canonical 
ensemble, see Eq. (124(1 . It would be interesting to see 
if the dependence of gap on system size could be deter- 
mined analytically for the microcanonical ensemble using 
random walk arguments. 

The distribution of the difference ln(A) — ln(Cii) + 1 is 
plotted in Figs.EJand^l In the canonical case, Eqs. 11' ">D 
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FIG. 10: Distribution of [ln(A) - In(Cii) + 1]/(2L) 1/2 for the 
microcanonical case. Equations 122H and (I25H predict that 
ln(A) — Ih(Cil) + 1 should be identically zero in the thermo- 
dynamic limit. The data seems to be tending towards this for 
large L. 

and (|20|l predict a Gaussian distribution with zero mean 
and standard deviation \/2L for large L. Fig. |5| shows 
that this works very well for the full range of sizes stud- 
ied numerically. Equations l|22l) and l|25|) predict that 
ln(A) — ln(Cn,) + 1 should be identically zero in the mi- 
crocanonical ensemble for L — ■+ oo and hence its distri- 
bution should be a delta function at the origin. Indeed 
the distribution in Fig. 1101 is narrow and sharply peaked 
at zero with a width which deceases as L increases, con- 
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FIG. 11: Probability distributions of Gil and mfmf, for sys- 
tem sizes L — 32 and 128 for the canonical ensemble. 
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FIG. 12: Probability distributions of Cil and mfmf, for sys- 
tem sizes L — 32 and 128 for the microcanonical ensemble. 



sistent with these expectations. 

Finally we look at the distributions of C\l and m\m s L . 
The scaling variables are 2 A = — \h{C\l) / L 1 / 2 and 
d = — ln(mf m s L )/L 1 / 2 and relevant plots are shown in 
Figs. II II and 1121 for the canonical and microcanonical en- 
sembles respectively. We see that in the canonical case 
the overall distributions of I and d are different at large 
arguments, but they match very accurately at small val- 
ues of the argument leading to the same behavior for the 
averages [Cii]£ v and [mf m s L ]% v shown in Fig. EJ In the 
microcanonical case, the overall distributions are quite 
similar but the agreement at small values of the argu- 
ment is not as good as in the canonical case. This leads 
to a greater difference, shown in Fig. \7\ between the aver- 
ages [Cil]™ c and [mfm£]™ c than in the canonical case. 



To summarize this section, for the canonical ensem- 
ble, the end-to-end correlation function [Ci/Jav falls off 
at criticality with a power of L, as predicted analytically, 
Eq. J22J>- However, for the microcanonical ensemble it 
falls off much faster, as a stretched exponential func- 
tion of distance. The average gap, [A] av falls off with 
a stretched exponential form at criticality in both en- 
sembles, with probably the same dependence on L. 

V. DISCUSSION 

In this paper we have looked numerically at the fi- 
nite size dependence of various quantities for the ran- 
dom transverse field Ising chain (RTFIC) at criticality, 
for both the canonical and microcanonical ensembles of 
disorder. For quantities that span the system, m s and 
C\l, finite size scaling appears, at first glance, to indi- 
cate different correlation length exponents for the two 
ensembles: v = 2 for canonical and v = 1 for the mi- 
crocanonical. However, in contrast to Igloi and Ricger 7 , 
we conclude that the correct interpretation is that the 
true scaling exponent is the same for the two ensembles, 
v = 2, but that the amplitude for the leading v = 2 piece 
is zero, in the microcanonical ensemble, for quantities 
that span the system and are thus sensitive to the micro- 
canonical constraint. Our reasons for this are two-fold: 

1. For the quantity we calculated that does not span 
the system, the "bulk" magnetization m, the same 
correlation length exponent v = 2 was found for 
both ensembles, indicating that there is a correla- 
tion length in the microcanonical ensemble which 
diverges with the larger exponent v — 2. Hence, if 
this is not seen for quantities which span the sys- 
tem, the explanation must be that the amplitude 
is zero, not that this larger length scale does not 
exist. 

2. The analytical expressions for the distribution of 
the surface magnetization m s , first obtained by 
Fisher 1 -, show that the scaling variable is 5L 1 / 2 
(canonical) and 5L 1 / 2 (microcanonical) demon- 
strating that the true correlation length exponent 
is v = 2 in both cases. 

Our interpretation of the data implies that the 
inequalitj^ v > 2/d is satisfied (as an equality) for the 
RTFIC, in contrast to the conclusion of Pazmandi et al. 3 . 

We have also looked at the energy gap between the 
ground state and first excited state at criticality. For 
both the canonical and microcanonical ensembles, a 
stretched exponential decay describes the data. For the 
canonical case, the exponent fi (the power of L in the 
exponential) is exactly 1/3, and our numerical data are 
consistent with fi = 1/3 for the microcanonical case too. 

The present model, the RTFIC, is integrable and thus 
relatively simple. In particular, the existence of the sim- 
ple analytical expression for the surface magnetization 
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7ji s , Eq.lJSJ, is surely related to the integrable nature of 
the model. Furthermore, the microcanonical constraint 
II hi — II Ji enters directly in this expression. Thus one 
can plausibly see how the constraint might affect quanti- 
ties which span the system and cause amplitudes for these 
quantities to vanish. However, for non-integrable mod- 
els, including models in higher dimensions, one would not 
expect the microcanonical constraint to enter in a direct 
way even for quantities which span the whole system. 
Thus it seems unlikely to us that there would be even 
an apparent difference in the critical behavior of non- 
integrable models in the two ensembles. We also note 
that the microcanonical and canonical ensembles of dis- 
order have been investigated for finite-T transitions in 
random systems by Aharony et al^ They find no dif- 
ference asymptotically between the critical behavior and 
finite-size effects of the canonical and microcanonical en- 
sembles (which they term grand canonical and canonical 
respectively) . 
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APPENDIX A 



where 



1 In hi — In Ji 
si = ~r~ 



L var(ln h) + var(ln J) 



(A2) 



Thus S is a sum of L uncorrelated random numbers 
i = 1, ...L with mean 



[£i]av — S/L 



and variance 



L z var(ln h) + var(ln J) 



(A3) 



(A4) 



Using the central limit theorem and the definition of £ in 
Eq. J2J we find that, in a canonical realization with given 
S the probability, Pg(S), of obtaining the precise value S 
is 



1/2 



P s( S )=[-) e-'V-V , (A5) 

for L — > oo. 

Now let P C (A, 5) and P mc {A, 5) be the probability dis- 
tributions of some observable A in the canonical and mi- 
crocanonical ensembles respectively. The two are related 
by 



P C (A,6) 



P mc {A,S)P s (6)d5. (A6) 



Measurements made in the two ensembles are in fact 
related to each other by a simple transformation at large 
L. To see this note that 



(Al) 



Correspondingly, expectation values in the two ensembles 
are related by 

/>oo 

M= [A]Z(S)P- s (5)dS. (A7) 
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